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2 Abstract. 

3 To investigate the internal structure of the magnetopause with spacecraft 

4 data, it is crucial to be able to determine its normal direction and to con- 

5 vert the measured time series into spatial profiles. We propose here a new 

6 single-spacecraft method, called the BV method, to reach these two objec- 

7 fives. Its name indicates that the method uses a combination of the mag- 

8 netic field (B) and velocity (V) data. The method is tested on simulation and 

9 Cluster data, and a short overview of the possible products is given. We dis- 

10 cuss its assumptions and show that it can bring a valuable improvement with 

11 respect to previous methods. 
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1. Introduction 

12 The Earth magnetopause is the outer boundary of the terrestrial magnetosphere. Cut- 
is side of this boundary, the magnetosheath plasma is the shocked solar wind plasma, i.e. 

14 cold and dense, with a magnetic field direction essentially determined by the solar wind 

15 one. Inside of it, the magnetospheric plasma is comparatively hot and tenuous, with 

16 a magnetic field direction essentially determined by the planetary one. Experimentally, 

17 investigating the magnetopause structure by spacecraft measurements is made difficult 

18 by the fact that the boundary is not steady: it can be shaken by the variations of the 

19 solar wind pressure, and perturbed by different kinds of waves, incident body waves as 

20 well as surface waves. It can also be locally and temporarily the place of different surface 

21 instabilities, implying or not magnetic reconnection, such as Kelvin-Helmholtz, Rayleigh 

22 Taylor or tearing instabilities {Hasegawa et al, 2012). 

23 Two informations are crucial to investigate the magnetopause nature: 1) accurately 

24 determine the direction of its normal with respect to the magnetic field (in a strictly 

25 stationary configuration, having the normal magnetic field i?„ null or not have quite 

26 different consequences on the physical nature of the layer, even if the non null -B„ is 

27 small) and 2) determine an approximate spatial coordinate along the normal, to be able 

28 to draw the spatial profiles of the different relevant parameters, in the boundary frame, 

29 i.e independently of the velocity at which these profiles are traversed by the spacecraft. 

30 Several methods have been developed for both of these two purposes. 

31 To study the large scale shape of the boundary, its motion and its orientation, multi- 

32 spacecraft methods have been developed, particularly for the ESA Cluster mission 
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33 {Paschmann and Daly, 1998 and 2008). These methods are essentially based on tim- 

34 ing differences between spacecraft and all rely on strong assumptions on the boundary: 

35 its form (plane or slightly curved at the scale of the spacecraft tetrahedron), its station- 

36 arity (constant profile and width, hereafter CTA for "Constant Thickness Approach", 

37 {Haaland et al, 2004), or its velocity with respect to the spacecraft (hereafter CVA for 

38 "Constant Velocity Approach", Russell et al 1983). Others are single spacecraft: they also 

39 rely on assumptions on the boundary properties such as planarity and stationarity, but 

40 they use in addition theoretical knowledge on the measured physical quantities, such as 

41 conservation laws. When using in particular the magnetic field data, the MVAB method 

42 (Minimum Variance Analysis on the magnetic field B, Sonnerup and Scheible, 1998) takes 

43 advantage that div(B) = 0, which draws -B„ = est in the 1-D case. Its variant MVABC 

44 (C for corrected) adds the constraint Bn = 0, using the additional information that the 

45 magnetopause normal component -B„ is generally close to zero, if not strictly zero. This 

46 allows to handle cases when two components are nearly constant and not a single one 

47 [i.e. when two eigenvalues of the variance-covariance matrix are small). When the mag- 

48 netopause can be supposed 1-D and stationary but when its thickness is small, making 

49 the kinetic effects non negligible with respect to the MHD ones, the experimental profiles 

50 have to be compared with the kinetic models of the tangential layers that can be found in 

51 the literature (see de Keyser and Roth, 1998, for a review of the first models of this kind, 

52 and Belmont et al., 2012, for the most recent one). The experimental method developed 

53 in this paper should enable to perform such comparisons. When the magnetopause layer 

54 cannot be supposed 1-D, other methods are needed. Some have been developed to recon- 

55 struct the magnetopause structure, supposing it is 2-D and stationary, and that it respects 
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56 MHD equations: these are the Grad-Shafranov reconstruction methods (see Hasegawa et 

57 al, 2004, for long-duration reconstruction). A review and discussion of short- and long- 

58 duration methods is made in De Keyser (2006). Experimentally, it is often difficult to 

59 decide whether the 1-D or the stationary hypothesis has to be questioned first. Future 

60 comparison between the results of the reconstruction methods and those of the method 

61 proposed in this paper should be interesting in this respect. 

62 To find an approximate normal coordinate allowing to investigate the internal struc- 

63 ture of the layer and to determine profiles across it, other methods have been developed 

64 independently, introducing the notion of "transition parameter" (Lockwood and Hapgood, 

65 1997). These methods can be used with single-spacecraft data. They also rely on assumed 

66 magnetopause properties, and they have been based hitherto on the variations in density 

67 and temperature of the electron population. This of course limits the temporal resolution 

68 of the method -and consequently its spatial one- to the electron experiment resolution. 

69 We propose here a new single spacecraft method, referred hereafter as "BV" to show 

70 that the magnetic field and the fiow velocity data are used simultaneously, to analyze 

71 magnetopause-like interfaces. It combines the two previous types in such a way that it 

72 allows to determine in the same operation the magnetopause normal with an improved 

73 accuracy and a transition parameter with an improved time resolution and expectingly 

74 closer to a real spatial coordinate. Fitting the magnetic field hodogram with a prescribed 

75 form, which is here an elliptical arc, allows to determine the normal direction with a fairly 

76 good accuracy. In addition, the angle a characterizing the position on the elliptical arc 

77 provides a reliable transition parameter inside the current layer, which can be viewed 

78 as a proxy for a normalized coordinate in the normal direction. On the other hand, as 
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79 soon as the normal direction is known, the velocity measurements give a non-normalized 

80 normal coordinate, which is just the integral of the normal flow velocity «„. It can give, 

81 in particular, a fairly good estimate of the physical width of the layer whenever the mea- 

82 sured velocity should be in most cases dominated by the motion of the boundary. Using 

83 simultaneously the magnetic and velocity measurements just consists in imposing that the 

84 normal coordinate determined by the only velocity measurements is proportional to the 

85 transition parameter coming from the only magnetic measurements. Since the integral 

86 of Un is very sensitive to the normal direction, this enables to improve the determina- 

87 tion of this direction with respect to the only magnetic one, while the time resolution 

88 of y{t) remains approximately the magnetic one, which is much better than the velocity 

89 resolution. 

90 Section 2 presents the principles of the BV technique, and section 3 the different valida- 

91 tion tests performed. The method allows to draw spatial profiles of any physical parameter 

92 across the magnetopause boundary. Examples of such profiles are presented in section 4, 

93 before discussing the interest and the limitations of the BV method and concluding in 

94 section 5. 

2. Principles of the method 

95 As the previous equivalent methods, the basic assumption of the BV technique is that, 

96 apart from oscillating perturbations, the boundary is sufficiently one dimensional and 

97 stationary at the scale of the spacecraft crossing. To explain the principles of this method, 

98 we use here a set of Cluster data on March, 3rd, 2008, when Cluster C3 encounters the 

99 magnetopause around 23:16, as it can be seen on Fig. 1 from the transition in the 

100 energy composition of the plasma, the density gradient and the rotation of magnetic field 
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101 observed. The method uses principally the magnetic field data {Balogh et ai, 1997). In 

102 subsection 2.1 we describe how we obtain an initial guess with only magnetic field data. 

103 Subsection 2.2 then explains the BV method itself, which combines magnetic field and 

104 ion velocity data. 

2.1. Initialization with the only magnetic field data 

105 In order to correctly initialize the minimization process of the complete BV method, 

106 involving magnetic field and ion velocity data, it is necessary to perform first an initializa- 

107 tion stage, which provides an approximated frame and a first elliptical fit. This stage uses 

108 only the magnetic field data. It is done itself in several steps. The first step consists in 

109 finding a first approximation of the normal direction via a MVABC technique {Sonnerup 
no and Scheible, 1998). Fig. (2) shows the tangential hodogram derived by this method. In 

111 this example as in many other observations {Panov et al, 2011), we observe a C-shaped 

112 hodogram, which can be fitted by an elliptical model. Although the general concept of 

113 the BV method is valid for any ID layer, its present implementation is conceived for such 

114 kind of hodograms. Further generalization to more complicated hodograms (in particular 

115 for the S-shaped hodograms described in Panov et al, 2011) is of course always possible. 

116 The second step consists in selecting the "magnetic ramp", i.e. the interval of data where 

117 the gradient of Bl is located. We then further select the data points by choosing only 

118 a sample of "representative points" among them. This step has a double purpose: elim- 

119 inate the perturbations that can be considered as "noise", and make the different parts 

120 of the crossing equally represented in the statistics, even if the spacecraft does not spend 

121 the same time in these different parts. First we roughly eliminate the perturbations by 

122 discarding all points too far from the mean trajectory of the hodogram, and we represent 
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each too close packet of points by only one single point. An elliptic fit and a new reference 
frame are derived from these points, using a Powell algorithm. The points selected in this 
way and the correspondent fit in the new frame are shown in Fig. (3). Then, the second 
goal is achieved by keeping a constant number of points in each a slice, which corresponds 
to the hypothesis (to be justified in next section) that a varies linearly with y. A new 
elliptic fit and approximated frame are then obtained, which provides a fine initialization 
for the BV method itself. 

2.2. Simultaneous use of magnetic and velocity data 

The above stage has given an initial guess for the BV method regarding 1) the normal 
direction, and 2) the parameters describing the elliptic hodogram. The main part of the 
method then consists in using the temporal information B(t), together with the velocity 
measurements from the Hot ion analyser experiment {Reme et al, 1997). Going back to the 
totality of the B data points, one minimizes the distance between them and the elliptical 
model B(|/), the function y{t) being the integral of the normal velocity Un- We therefore 
assume that this velocity is dominated by the layer velocity, i.e. that the normal velocity 
in the layer frame is zero or negligible. The minimization is done with respect to the three 
angles that characterize the rotation of the ellipse proper frame and to the parameters of 
the elliptic hodogram, initialized previously, using the same Powell algorithm as above. 
The distance to be minimized is: 



Z2 yi^dx — BmxY + {Bdy — B„iyY + {Bdz — B^zY (1) 
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Where B^ represents the data points and B^ represents the model. This model is given 
by: 

Bmx = B^ocosa (2) 

Bmy = ByQ (3) 

Bmz = B^osina (4) 

with: 

a = ttl + (a2 - ttl) y/Vmax, (5) 

y being the position deduced from the normal velocity integral. The magnetic field data 
and velocity data are obtained from prepared data by a rotation of M{6,(f),x)- The 
parameters of the fit are 9, 0, x, B^o, ByQ, B^o, «!, and a2- 

This final stage provides all the needed outputs: the normal direction, the spatial 
position y{t) along this normal (measured directly in physical units, providing in particular 
the layer thickness in km), and the fit of magnetic field, as illustrated, for example, on 
Fig. (4). Here the computed magnetopause thickness is 1800 km and the linear Pearson 
correlation coefficients of the fit of B^ and Bz are 0.99 and 0.95. The spatial position y 
is then extrapolated linearly outside the boundary, in order to plot approximated profiles 
of any plasma parameter on scales larger than the ramp region if necessary. 

3. Validations of the method 

Having presented how the BV method works in the previous section, we will now explain 
what led us to this way of proceeding and what are the different validation tests we have 
performed. We will discuss first the validity and the limitations of the hypotheses done, 
and discuss afterward the consistency of the obtained results. We used three different 
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tools to develop and validate the method: - a simple code to generate artificial magnetic 
field data, - a hybrid simulation of an asymmetric reconnection layer {Aunai et al, 2013b), 
- and real data from the Cluster mission, especially a 2008 low latitude crossings list 
compiled by N. Cornilleau-Wehrlin. 

3.1. Hypotheses: elliptical shape and linear angular velocity 

The first new assumption of the method, with respect to previous single spacecraft data 
analysis methods, is the elliptic shape of the tangential magnetic field hodogram, the 
simplest model geometry to describe C-shaped hodograms. This elliptical shape is indeed 
consistent with a simple generalization of the circular model B(?/) proposed by {Panov et 
al, 2011): 

Blo 

Em 1 



tanh(|//L) (6) 

(7) 



Bmo cosh{y/L) 

These formulas imply in particular that B\/ B\q + B\j/ B\,iq = 1, which can be a test of 

the elliptical shape. 

The efficiency of the method can be tested first on a numerical simulation of reconnec- 
tion {Aunai et al, 2013b), far from the X point. Its applicability is not obvious in this case, 
since, before the development of the reconnection pattern, the initial condition is purely 
tangential, without any rotation. Nevertheless, the Hall effect creates a self-consistent 
out-of-plane magnetic component during the reconnection process, which, in the consid- 
ered asymmetric configuration (asymmetric in density and temperature and coplanar and 
antisymmetric in magnetic field), results in a C-shaped hodogram if looked between the 
separatrices. Fig. (5) shows the magnetic field in the interval that corresponds to the 
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gradient of Bl. The error is here less than 2 percent. We have checked that this good 
accuracy is kept as long as the crossing considered is not too close to the X-point, which 
is generally the case for crossings of reconnected magnetopause or to the limits of the 
simulation. 

Concerning the analytical form of a{y), we also checked the validity of the linear hy- 
pothesis in the same simulation study. Fig. (6) shows how a varies as a function of 
the normal coordinate ys of the simulation . We observe that, apart from weak periodic 
variations, the linear form is well satisfied. It is worth explaining that the weak periodic 
departures from the linear variations (which can be well described by the three of four first 
terms of a Fourier transform) can indeed be accounted for in the minimization procedure, 
but it would increase the number of free parameters and drastically affect the convergence 
of the minimization process. 

3.2. Consistency of the results and limitations 

Regarding the consistency of the results, the first test consists in running the first part 
of the method (identification of the ellipse and of its proper frame) on a magnetic field 
that is artificially generated with an elliptic hodogram. Such artificial data have thus 
been constructed with the same analytical formulas as those of the program, then turned 
on a random frame, and added with a random Gaussian noise centered on the signal, 
with a relative amplitude up to 50%. The result is that the method always allows to find 
the good initial normal direction with at least 5 significant numbers, as well as the right 
ellipse parameters, whenever the noise does not exceed 30%. 

The second test consists in using the above numerical simulation {Aunai et al, 2013b) to 
mimic a real magnetopause crossing. In order to make the method work, we must modify 
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the simulation results in a way that makes it likely closer to most real magnetopause 
crossing: we multiply the tangential velocities by a large factor {^ 10). Thanks to this 
change, the tangential velocities get a much larger contrast than the normal ones, which 
is necessary for the program convergence. It must be noted that such a contrast of 
the tangential velocities does generally exist at the magnetopause, since the tangential 
velocity change is generally of the order of a few 100 km/s, while the normal one (in 
the spacecraft frame) is generally about ten times smaller and varies very little. In order 
to focus on the reconnection process freely of any KH instability, the simulation did not 
include such a velocity shear. Furthermore, the normal velocity of the virtual spacecraft 
considered with respect to the boundary, has to be chosen large enough with respect to the 
normal velocities in the boundary frame. This is also, as already mentioned, a reasonable 
hypothesis for a real magnetopause crossing. 

Under these assumptions, we get normals with an angular precision oscillating between 
and 5 degrees (with the corresponding errors on the shape of the tangential hodogram) 
and 0-5% errors on the y parameter (and derivative), which corresponds to the internal 
velocity and the approximations on a{y). 

The result is not changing as long as the virtual spacecraft crosses the simulation far 
enough (several di) from the X point, where the 2D effects are not dominant. In these 
cases, the precision of the MVABC method is of the same order, (slightly better or worse, 
depending on the cases), because -B„ is actually very close to 0. 

Regarding real Cluster data, the measurements show more perturbations, but the varia- 
tions of the field value around the mean ellipse are still around 5 percent for most C-shaped 
hodograms. A good test for the elliptical shape is to plot B^^B^), that should be linear 
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210 for a tangential ellipse. Fig. (7) shows this plot for two magnetopause crossings on 

211 03/03/2008 and 04/01/2008. It shows that the elliptical shape is a good approximation. 

212 It is clear, from the the tests on the numerical simulation, that the BV method has 

213 limitations related to the necessary contrast between the normal and tangential component 

214 profiles. When applied to real Cluster data, these limitations may have, in some occasions, 

215 consequences on the results obtained. We will discuss these limitations in the conclusion 

216 section. It is to be noted however that these limitations are based on assumptions which 

217 are different -and generally weaker- than those of the other single-spacecraft methods such 
2ia as MVAB or MVABC. 

219 We will present a detailed study on a case {Dorville et al, 2013), where the BV method 

220 leads to a better understanding and more precise results than MVAB(C). When all the 

221 methods are confidently applicable, the results seem to be consistent with each other and 

222 with the theoretical knowledge. We show on Fig. (8)" oeil" a reproduction of a figure from 

223 {Haaland et al, 2004) corresponding to a benchmark case where different methods have 

224 been used. The center of the figure is the mean MVABC normal and other single and 

225 multi-spacecraft methods are represented in a polar plot in the plane perpendicular to this 

226 normal. The result of the BV method on CI spacecraft is indicated by a star. The figure 

227 shows that, if the result is different from other methods, it is inside the dispersion range of 

228 the points. The thickness of the layer always stands between a few hundreds of kilometers 

229 and a few thousands, which is consistent with literature, the tangential velocities being 

230 generally one order of magnitude larger than the normal one (in the spacecraft frame). 

231 The normal magnetic fields always stand between and 20 nT, the non null values being 
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232 reliable and quantitative indications of a connected boundary, which could hardly be 

233 obtained previously. 

4. Products of the method 

234 As explained above, the first main direct product of the method is an accurate deter- 

235 mination of the direction normal to the boundary, leading to reliable values of the small 

236 components -B„ and «„ of the magnetic field and the flow velocity across the boundary. 

237 The second direct product is the determination of a spatial coordinate y{t) allowing to 

238 draw any plasma parameter profile against the spatial position y from their temporal mea- 

239 surement. The magnetopause layer thickness is also an interesting by-product deriving 

240 directly from the two preceding ones. 

241 Examples of y profiles are presented in Fig. (9) for the crossing of 03/03/08. Here we 

242 see the characteristic jump of density at the magnetopause, but no temperature jump, 

243 the pressure evolving like the density. For the different crossings that we investigated, 

244 we could often observe clear differences concerning the locations of the particle gradients 

245 with respect to the magnetic field rotation. In a companion paper, we will present an 

246 interesting case study where the BV method can bring new information about the nature 

247 of the magnetopause. 

248 In Fig. (10) the normal electric field obtained with the EFW experiment {Gustafsson et 

249 a/,(2001)) and the tangential components are shown for the same 03/03/2008 crossing. We 

250 see that the maximum variance is on the normal electric field, as expected by theory, and 

251 quite constant tangential electric fields, which confirms that the normal direction found 

252 is a good one. Fig. (11) shows the profiles of magnetic field spectral power density ob- 

253 tained with the STAFF experiment {Cornilleau-Wehrlin et al, 2003)for different frequency 
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254 ranges. One can observe that the source of waves hes in the magnetosheath and that the 

255 depth of penetration depends on the frequency, the lowest frequencies penetrating deeper 

256 toward the magnetospheric side. 

257 This abihty to get spatial profiles of all the quantities in the boundary is a key to a 

258 better understanding of the physical nature of the magnetopause. 

5. Discussion and conclusion 

259 We have presented the new BV method to analyze the structure of the magnetopause 

260 boundary layer, using spacecraft data. It combines the magnetic field and velocity mea- 

261 surements of one single spacecraft and permits to find the normal direction and a good 

262 resolution on a spatial coordinate to resolve small scale variations inside the layer. Using 

263 it, we are able to study the internal structure of the layer, for any of the physical quan- 

264 titles measured on board. The method works on simulation and generated data, and its 

265 assumptions can be verified on Cluster crossings. 

266 It is worth observing the conditions of validity of the BV method are not the same 

267 as the other single spacecraft methods such as MVAB, and that they are in general less 

268 restrictive. In MVAB, one needs to discriminate Bn and Bm, which fails systematically in 

269 structures as shocks, and often at the magnetopause since this one is often quasi-coplanar. 

270 MVABC has the same condition of validity, with the additional problem that it cannot 

271 be used for determining i?„ since this component is supposed null. In the BV method, 

272 one needs to discriminate the two couples of data sets: (-Bat, Vn) and {Bm, Vm)- This is 

273 clearly a weaker condition since, even if Bn and Bm are nearly constant, the differences 

274 between Vn and Vm, (profiles and/or orders of magnitude) are generally sufficient to 
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275 guarantee a correct operation. The difficulties can only arise when not only Bn and Bm 

276 are indistinguishable (mean jump much smaller than noise), but also V^ and Vm- 

277 Contrary to the multi-spacecraft timing methods, the BV method can also handle cases 

278 when the boundary is shaken with a non trivial normal velocity evolution (which seems 

279 frequent). When this evolution is non negligible between two spacecraft crossings, the 

280 timing methods obviously fail. 

281 The BV method however brings a new limitation: although one works essentially with 

282 magnetic field data, a sufficiently long crossing is needed (at least three or four velocity 

283 measurement points inside the crossing)to make efficient the contribution of the velocity 

284 data. We are therefore not able to analyze as many crossings as the other methods. 

285 With the proposed method, the structure of the magnetopause should be now open to 

286 more detailed investigations. Some examples of spatial profiles have been given in section 

287 4. The method is used in a companion article, for an atypical magnetopause case study 

288 giving new insight on this structure. 
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Figure 1. Density, energy spectrogram and magnetic field observed by Cluster CI around 
23hl6 on 03/03/2008. The jump of density, change in plasma energy composition and rotation 
of magnetic field show that the satellite is crossing the magnetopause. 
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Figure 2. Hodogram of the magnetic field in the tangential MVABC plane for the Cluster C3 
magnetopause crossing of 03/03/2008. 
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Figure 3. Initialization fit of tlie liodogram in initialization frame for the Cluster C3 mag- 
netopause crossing of 03/03/2008. The data is in black, the selected points in yellow and the 
initialization fit in red. 
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Figure 4. Fit (dashed) of the three components of the magnetic field (solid lines) along the 

normal coordinate for the Cluster C3 magnetopause crossing of 03/03/2008. 
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Figure 5. Fit of tlie liodogram of tlie magnetic field in the simulation for a virtual satellite 
crossing far from X point in the numerical simulation. The data is in black, the selected points 
in yellow and the fit in red. 
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Figure 6. Angular position on the ellipse along the normal direction in the simulation for a 
virtual satellite crossing far from X point. 
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Figure 7. Hodograin of the squared components of magnetic field in the LM plane for two 
Cluster C3 magnetopause crossings of April, 1st 2008 and March, 3rd, 2008. 




Figure 8. Several single and multi-spacecraft methods normal direction positions in the plane 
perpendicular to the MVABC mean normal for a benchmark case from {Haaland et al, 2004). The 
star represents the result of BV method. The distance from the origin in this plane corresponds 
to the sin 6, where 6 is the angle from the mean normal direction. 
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Figure 9. Evolution witli normal position of density, pression and temperature measured by 
Cluster C3 for the 03/03/2008 crossing. 
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Figure 10. Evolution with normal position of the Electric field measured by Cluster C3 for 
the 03/03/2008 crossing. 
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Figure 11. Evolution witli normal position of the spectral power density for magnetic field 
measured by Cluster C3. The total spectral power density is the thick line, as the daHz hHz and 
kHz frequency ranges are also represented on the same scale. 
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